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We prove that sign problems in the traditional approach to some lattice Yukawa models can be completely 
solved when fermions are formulated using fermion bags and bosons are formulated in the worldline repre- 
sentation. We prove this within the context of two examples of three dimensional models, symmetric under 
Ul,(1~) X Ur(1) X (Parity) transformations, one involving staggered fermions and the other involving Wilson 
fermions. We argue that these models have interesting quantum phase transitions that can now be studied using 
Monte Carlo methods. 
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I. INTRODUCTION 

Feynman path integrals can be used to map a quantum sta- 
tistical mechanics partition function into a classical statistical 
mechanics partition function with one caveat : the Boltzmann 
weight of the classical partition function may be negative or 
even complex. When this occurs the mapping is said to suf- 
fer from a sign problem since the mapping is not useful for 
Monte Carlo methods. However, fortunately the mapping is 
not unique and it may be possible to find a different map- 
ping in which the Boltzmann weights of the classical partition 
function are indeed positive. If these weights are calculable 
with polynomial effort as the system size grows, the mapping 
is said to be free from a sign problem and one may be able 
to construct a Monte Carlo algorithm to sample the classical 
configuration space and solve the problem. In many interest- 
ing cases a mapping without a sign problem has eluded re- 
searchers and discovering the correct mapping is defined as a 
solution to the sign problem. 

Sign problems are well known obstacles to solving many 
quantum field theories from first principles. This is partic- 
ularly true when the microscopic degrees of freedom con- 
tain fermions. Famous examples of physical systems where 
sign problems have hindered progress are a finite density of 
strongly interacting matter |Q}[3) and a finite density of elec- 
trons with Coulomb repulsion |4|. Many other model field 
theories also suffer from sign problems. Examples of these 
include four-fermion field theories and interacting boson- 
fermion (Yukawa) models. These often have interesting low 
energy physics and are used as effective field theories [5]|6). In 
three dimensions they can contain interesting quantum phase 
transitions that have remained unexplored due to sign prob- 
lems. 

While some important sign problems may be unsolvable 
0, it has become clear that a variety of sign problems can in- 
deed be solved by finding the right representation of the par- 
tition function. Novel solutions to sign problems have been 
found in both purely bosonic models with a complex action 
ll8 ffl"3l and in purely fermionic models lfl4l[T5l . Here, for the 
first time we show that sign problems in models containing 
fermions and bosons as dynamical fields interacting with each 
other (which we call Yukawa models) can also be solved. We 
illustrate this using two examples of lattice models in three 
dimensions. 



In section 2 we review the sign problems that haunt the tra- 
ditional approach to two lattice Yukawa models one with stag- 
gered fermions and one with Wilson fermions. In section 3 we 
show how the fermion-bag approach when combined with the 
worldline formulation of bosonic degrees of freedom solves 
these sign problems. In section 4 we present our conclusions. 



H. SIGN PROBLEMS 

In this section we review two examples of sign problems in 
lattice Yukawa models, one with staggered fermions and the 
other with Wilson fermions. For convenience our models are 
defined in three space-time dimensions, but the discussion is 
applicable in higher dimensions with minor modifications. We 
begin with the example with staggered fermions. The action 
of the model is given by 

S s = J2^x(D S [9]) xy %+S b [e] (1) 

where i[> x , tf'x are two Grassmann valued fields on the lattice 
site x = (xi, X2, £3) of a cubic lattice with V sites. While 
there are many choices for the bosonic action, for simplicity 
in this work we choose it to be the classical XY model, 

S b [e} = -j3j2 c °s(6 x -0 y ). (2) 

where the bosonic field 8 X is a phase. Here (xy) refers to near- 
est neighbor sites. The matrix D 3 [9) is the V x V staggered 
Dirac operator whose matrix elements are given by 

(D s [d]) xy = -g e ie * e *5 x , y + (D s °) xy , (3a) 

(I>* )*V = $>a,xV£,, (3b) 
a 

where the fluctuating mass term depends on the bosonic field. 
The index a = 1,2,3 represents the three directions, f] a ,x 
are the staggered fermion phase factors (r/i = 1,772 = 
{~l) Xl , m = (-l) Xl+X2 ), e x = (_l)*i+*»+*3 i s t he site 
parity and 

^xy — 2^ x >y+ & _ $x+&,y)- (4) 
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FIG. 1. A schematic phase diagram of Ul(1) x Ur(1) x Z2 sym- 
metric lattice Yukawa models discussed in the text. The symmetric 
phase contains massless fermions while the broken phase contains 
massless bosons. The critical point on the g — line is the XY 
critical point. For small f3 the model reduces to a Gross-Neveu (GN) 
model. 



The parameters f3 and g control the physics of the model. It is 
easy to verify that the action in Eq. ([l} is invariant under the 
following Ul(1) x Ur(1) chiral transformations : 

^ e i( Mi+^)/2+iMi-e*)/2 ^ (5a) 
^ e -ie L {i-^)/2-ie R (i+e x )/2 ^ (5b) 

X -> 6 X + (6 R - 9 L )e x (5c) 

In addition, the action is also invariant under the Z 2 (parity) 
transformations: 



i> x -> -tp- x , tpx iP- 



— ► 



+ TT (6) 



These symmetries play an important role in governing the 
long distance physics of the model. The model possesses at 
least two phases, a symmetric phase with massless fermions 
and massive bosons and a broken phase with massless bosons 
and a massive fermion. Since a fermion mass term breaks 
parity in three dimensions, it is natural to expect massive 
fermions in the broken phase. Both the phases can be accessed 
by tuning j3 and g. A schematic phase diagram, expected from 
general arguments, is shown in Fig.[T] The critical point on the 
g = axis is the well known 3d-XY critical point. If it exists, 
a second order quantum critical line separating the two phases 
starting from the XY point must be governed by a different 
critical point. Based on symmetries we conjecture that it be- 
longs to the universality class of the model studied in lfl6l . 
Unfortunately, as we will discuss below, this lattice Yukawa 
model cannot be studied using traditional Monte Carlo meth- 
ods due to sign problems. 

Next we consider a Yukawa model constructed with Wil- 
son fermions. We again restrict ourselves to three space-time 



dimensions for simplicity. The action of the model is given by 



where, unlike the staggered fermion case, the fields ip x and 
ip x are each four component fields, written in terms of four 
two-component left and right fields as 



(8) 



We label the eight Grassmann fields on each site as ijj g x , ip g x 
where c = L. R and s = 1,2. The Dirac operator is given by 



(D w °) xy -gle~ ie 
gle ie *6 xy (W) 



(9) 



where a are the Pauli matrices and I is the 2x2 identity ma- 
trix. The 2V x 2V matrix D w ° is the Wilson-Dirac operator 
defined by 



(D Wo ) xy =IM xy +<j-V x 



(10) 



where V was defined in Eq. Q and M is the V x V Wilson 
mass matrix defined by 



M X y — Xt y , ^ 



y+a 



The action ^ is also invariant under the following Ul(1) 
Ur(1) transformations: 



% + {e R -e L ). 



(12a) 

(12b) 
(12c) 



However, due to the presence of the Wilson mass term M xy , 
the action is not invariant under the Z2 (parity) transforma- 
tions 

x^-x, ip x -ip_ x , ip x ^4>- x , 9 X ->• 0- x + ir. (13) 

Parity transformations are restored in the long distance 
physics, when the hopping parameter n is tuned to a critical 
value where massless fermions emerge. For g = 0, one finds 
k c = 1 /3. In general k c is a function of the couplings g and (3. 
Once k is tuned to this value the phase diagram of the model 
as a function of j3 and g is very similar to Fig. [T] except that 
the critical behavior on the quantum critical line could belong 
to a different universality class as compared to the staggered 
fermion case. The difference is essentially due to the number 
of fermions that become massless on the critical line. 

In order to study the two Yukawa models discussed above 
one begins with the partition function 



Z t = / [# dip] 



(14) 
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where i = s, w and tries to construct a Monte Carlo technique 
to compute expectation values of appropriate quantities. The 
traditional method is to integrate the fermions completely and 
rewrite 



[dO] e 



Bet{D l [9}). 



(15) 



If the determinant of the Dirac operator D l [6] was non- 
negative then one could have devised a Monte Carlo method 
to sample the [6] configurations. Unfortunately this is not the 
case for both the models considered here. With staggered 
fermions there is no symmetry that can be used to even ar- 
gue that Dct(D s [0]) is real. The determinant can in fact be 
complex. On the other hand with Wilson fermions one can 
prove that 



C{D W [9])C = (D w [0])*, where C 



— 1(72 

ia 2 



(16) 

which means the determinant of D w [9] is real. Unfortunately, 
it could still be negative. Thus, both Yukawa models discussed 
above suffer from sign problems in the traditional method. As 
far as we know, these sign problems have remained unsolved 
until now. While these sign problems clearly arise for reasons 
completely different from the sign problem in finite density 
QCD J3], they seem equally difficult. 



III. FERMION-BAG WORLDLINE APPROACH 

We now prove that the sign problems disappear in the above 
two models when the fermions are formulated in the fermion- 
bag approach and bosons are formulated in the worldline rep- 
resentation. Since the details of the proof are slightly different 
in each case, we will discuss them separately. Let us first con- 
sider the partition function of the Yukawa model with stag- 
gered fermions. As explained in |[T6l[T7l we expand the inter- 
action in powers of the coupling. We begin by noting that 



= l+ff(lM>«)e* 



(17) 



due to the Grassmann nature of ipxi^x- We can represent the 
two terms on the right through the discrete variables n x — 
0, 1. The first term refers to n x = or no interaction while the 
second term represents n x = 1 and indicates the presence of 
an interaction vertex ip x ip x (or a monomer). Every monomer 
is also associated with the term e i£x9a! . Using this idea it is 
possible to write 

[n] 

x J [dO] e~ Sb[e] e fe "i*"i...e te "* fl "J } (18) 

where [n] represents a configuration of monomers, j refers to 
the total number of monomers and 



(19) 



is the free staggered fermion action. . For every configuration 
fnl we have labeled the sites where the monomers are located 



as { 
that 



j}. As explained in 1 17, 18 1 it is possible to show 



# #] e" s ° ^ x ^ xx ..4 x .^ Xj = Det(W[n}) > (20) 



where W[n] is a (V—j) x (V —j) matrix which has the same 
elements as D s ° except that the sites xi,..,Xj are dropped. 
The sites (V — j) form what we call the strong coupling 
fermion bag lfl6l [TSI . Interestingly, one can also rewrite the 
bosonic integral using the worldline representation as dis- 
cussed in lTTTl|T9ll . One finds that 



[dO] e 



-S h [0] e ie xi 9* 



[fc] x,a 



+ e x n x ^j (21) 

x a 

where the integer bond variables fc x , Q represent worldlines of 
charged particles, Ik x a (ft) is the modified Bessel function, [fc] 
represents a configuration of these worldlines. Note that the 
interaction between the fermions and bosons appear through 
the constrained delta function which essentially implies that 
every monomer either creates or destroys bosonic particles. 
Thus, using a combined fermion-bag-worldline representa- 
tion, the partition function of the staggered Yukawa model 
takes the form 

Z S = J2 .9 J Det(W>D H{Ik., a (P)} 



[n,k] 



| | fi( ^~^ (fcg,q k x — a , a) + £x n x^j 



(22) 



The Boltzmann weight of every [n, fc] configuration is non- 
negative and the sign problem is absent. 

Let us now turn to the lattice Yukawa model with Wil- 
son fermions. Unlike the staggered case, we now have four 



—L 



—L 



-tR 



iO, 



couplings gip hx ipi >x e w » , g^ 2 ,xW,x e > -9i>i,J>\,^ 

and — gtp 2 x^kx^*- We expand in each of these couplings 
and write 



r e ie x 



(23) 



for each of the four couplings. We then introduce four types 
of monomers ni t x) n 2,x> n 3,x, n 4,x = 0, 1 at each site x repre- 
senting these four couplings. For every monomer configura- 
tion [n] we label the sites where n\ >x , n2, x ,n>3,x an d n 4.x are 
non-zero, as Wi, Wj t , x%, xj 2 , y\ , . . . , y j3 , and z±, Zj A 
respectively. Thus, in the fermion-bag approach we can write 

z w = J2 g jl+j2+h+j A J [##] e~ s ° (-l) h+ji 



VI" 

I [d9] e~ Sh 



x,y 



"i...e ''...e "...e 



(24) 
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where 



is the free Wilson fermion action. Based on the fermion inte- 
gral it is easy to see that only those monomer configurations 
which satisfy j = ji + j 2 = 33 + 3 a contribute to the path 
integral. If we define the free fermion propagator 



Gs, x -s',y = J [d^d^j] e s °^ x ip s , ty 
then it is easy to prove that 

{G s ,x;s',yT = / e~ S ° (-1)^^. 



R 

«' ',V 



(26) 



(27) 



Using this result along with Wick's theorem it is possible to 
show that 



x ^2,. 1 ^^--V'2,^ 2 V'2\ 2 V'i; yi V'f, wl ---V'i, W3 V'i L yj3 



2,Zj 



|Det(G[n])| : 



(28) 



where G([n]) is a j x j matrix whose matrix elements 
are the free fermion propagators defined in Eq. ( [26] ) from 
(s = 1, {wi, ...ajjj}) and (s = 2, {xi, x j2 }) to (s 
1 ■ {t/i, -t/ja}) and (s = 2, {zi, z u }). Combining this re- 
sult with the worldline representation of the bosonic integral, 
(similar to Eq. p2]>), we see that 



Z W = J2 9 2J |Det(G[n])| 2 ]J{I K JP)} 



[n,k] 



]]^(t, tt - k x - a ,a) - ni iX - n 2 , x + n 3>x + n 4 . x ) 



(29) 

The above expansion of the partition function is again free 
of any sign problem since the Boltzmann weights are non- 
negative. 



IV. DISCUSSION 

In this work we have shown that the fermion-bag approach 
along with the worldline representation of bosonic degrees 
of freedom, allows us to solve sign problems in some lattice 
Yukawa models. While the solutions presented here depend 
on the details of the models considered, the idea is more gen- 
eral and can be applied to solve other sign problems, including 
those in non-relativistic field theories. For example, models 
that contain pairing interactions of the form gip^tpf-e can 
also be solved. Further, we can study models with discrete 
rather than continuous symmetries. For example, if we re- 
place Eq. pa} by 



(30) 



the C/l(1) x Ur(1) symmetry of the model is reduced to a 
Uf(l) x Z-2- While the conventional approach still suffers 
from a sign problem, the fermion bag approach does not. Sign 
problems can also be solved with Nf fermion flavors as long 
as all fermions couple to boson fields as in the one flavor case. 
This enhances the symmetry of the models by an SU(Nf) 
factor. 

The bosonic action can also be changed. For example in- 
stead of the XY model action, one can also choose the more 
standard form where 



s b W\ 



(31) 

In this case the phase e l6m in the couplings is replaced by a the 
complex field (j> x itself. Then the bosonic integrals in Eqs. ( 18 1 



and (24 1 can again be represented in a worldline represen- 



tation with positive Boltzmann weights. Indeed, by writing 
4>x = Px^ 6 " in the polar form one can show that 



= I [dp] e -£.C«*+ApJ) p Wl ..p W]Pzi ..p Z] 

X < II ih^AZPxPx+a)} 

I [k] x,a 

kx-a,a) + n x j > (32) 



x a 



where n x 



1 for x 



Wi, 



and n. 



-1 for x 



{D s [9]) xy = -g sin(6 x )5 x<y + (D s °) a 



Zi,..,Zj and n x = otherwise. Thus, the bosonic integral 
can again be represented without a sign problem. 

An important lesson from our work is that solutions to 
sign problems only emerge when along with pairing in the 
fermionic sector, the correct bosonic variables can also be 
identified. Sign problems in theories with multicomponent 
boson fields continue to remain a challenge, although our 
ideas may provide hints for finding a solution. In certain 
cases we believe that positivity may emerge in terms of quan- 
tities like fermionants which unfortunately can be exponen- 
tially difficult to compute [20|. The ability to convert a quan- 
tum partition function into a classical partition function should 
ultimately be dependent on the underlying physics and the de- 
grees of freedom that capture it optimally. 
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